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Abstract 

Hopping parameter expansions are convergent power series. Under general 
^< conditions they allow for the quantitative investigation of phase transition 

and critical behaviour. The critical information is encoded in the high order 
coefficients. Recently, 20th order computations have become feasible and 
used for a large class of lattice field models both in finite and infinite volume. 
They have been applied to quantum spin models and field theories at finite 
temperature. The models considered are subject to a global Z2 symmetry or 
to an even larger symmetry group such as O(N) with N > 2. 

In this paper we are concerned with the technical details of series com- 
putations to allow for a nontrivial vacuum expectation value < p(x) >^ 0, 
which is typical for models that break a global Z2 symmetry. Examples are 
scalar fields coupled to an external field, or manifestly gauge invariant effec- 
tive models of Higgs field condensates in the electroweak theory, even in the 
high temperature phase. A nonvanishing tadpole implies an enormous pro- 
liferation of graphs and limits the graphical series computation to the 10th 
order. To achieve the hopping parameter series to comparable order as in 
the Z2 symmetric case, the graphical expansion is replaced by an expansion 
into new algebraic objects called vertex structures. In this way the 18th order 
becomes feasible. 
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1 Introduction 



Linked cluster expansions (LCE) provide a way to construct hopping parameter or 
high temperature series of correlation functions for a large class of lattice models |], 
. In contrast to standard perturbation theory, the series are convergent power series 
(for a recent proof cf. || || and references therein). Under general conditions the 
domain of convergence extends to the phase boundary. The nature of the transition 
such as its location, the order and the critical exponents, are encoded in the high 
order behaviour of the expansion coefficients. 

In the recent past, much effort has been done to prolong the series for various 
field models on the lattice [Q-[|TUJ. Fast and efficient methods have been developed 
such that susceptibilities and renormalized coupling constants are available up to 
the order 20 in the hopping parameter. In particular, they have been generated 
for lattice field theories at finite temperature |J or quantum spin models. Partic- 
ular applications concern the universality classes of high temperature transitions 
for 0(N)-symmetric field theories |7j and Gross-Neveu models, related to the QCD 



chiral phase transition |l 1[ . More subtle questions concerning the applicability of 
dimensional reduction have been answered. The series have been generated and 
used in finite volume together with finite size scaling analysis to resolve the order of 
phase transitions Jl2|. 

The hopping parameter expansion amounts to a graphical representation of cor- 
relation functions. Appropriate graph classes have to be defined and generated, and 
their weights have to be computed. To make high order computation feasible, graphs 
are represented algebraically by incidence matrices. This representation needs to be 
one-to-one and fast ||. 

To achieve high orders, some properties of the model and of the lattice are taken 
into account in the course of the graph construction. Beyond others, it is used that 
the models considered are subject to a global Z 2 symmetry or a symmetry group 
including Z 2 as a subgroup. This implies that only even correlation functions are 
non-trivial. 

There are interesting models both in statistical mechanics and in quantum field 
theory where this global Z 2 symmetry is absent on physical grounds or by explicit 
breaking. We give two examples. 



• A scalar field theory on the lattice A, coupled to a nonvanishing external field, 

s(*) = -k HxMy) + Y,(^) 2 + K^) 2 ) 2 )+Hj2H^, (i) 

where the prime denotes nearest neighbour interaction, and H 7^ 0. Putting 
the hopping parameter k to its critical value k c (h, A, H = 0) and letting H — > 0, 
the critical exponent 5 is determined by the behaviour of < >. Knowledge 
of its series expansion in k allows for a high precision measurement of 5. 
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Effective models in which the fluctuation fields are confined in value to some 
subset of the real numbers. An example is provided by the gauge invariant 
SU(2) Higgs model. At finite temperature this model is of great importance 
for the physics of the early universe. The high temperature phase transition 
is highly non-perturbative in the region of realistic Higgs masses. On a finite 
temperature lattice A, integration of the compact gauge fields without impos- 
ing a gauge fixing generates an effective model that is described by a partition 



function of the form 13 



poo 

z= nw^w^t-w))- w 

The Higgs field condensates and their correlations are obtained as correlation 
functions of the field Although both the measure and the action V e ff are 
formally invariant under the global Z2 transformation £(x) — > — £(x), the field 
£ is confined to non-negative values, which implies that < >^ even in 
the high temperature phase. 

In this paper we generalize the linked cluster expansion to models in which the 
Z2 symmetry is explicitly broken or does not exist at all. There are mainly two 
generalizations to be included. First, correlation functions with an odd number 
of fields do no more vanish in general. In turn, vertices are allowed to have odd 
number of lines attached. This is easily incorporated into the known algorithm and 
requires only a slight generalization, although the number of graphs is considerably 
enhanced. Second, as the more subtle point, the vacuum expectation value of the 
field, the tadpole moment, does no more vanish in general, 

< p(x) > + 0. (3) 

This property requires the introduction of new graph classes with different topo- 
logical properties. For instance, one-particle irreducible (1PI) correlation functions, 
represented as a sum over 1PI graphs, play an important role in the course of the 
expansion. Once (Q) holds, however, 1PI refers only to external line channels. The 
necessary modifications will be given here. 

Furthermore, the combinatorical complexity is drastically enhanced. We expect 
that without further improvement the complexity to 10th order is already as much 
as the 20th order before. 

To achieve comparable order as in the Z 2 symmetric case (order 18-20), we 
replace where required the concept of graphs by the concept of vertex structures. 
Vertex structures have already been introduced in ||. They were used mainly to 
avoid a recomputation from scratch of LCE series if the coupling constants of the 
model are changed. Using them also avoids severe roundoff errors due to the finite 
precision of real coupling constants. 

Here, in addition we use vertex structures to avoid the graphical construction of 
the tadpoles and the combinatorical problems and computational limitations related 
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to them. In a way, we can view a vertex structures as an equivalence class in the set 
of graphs considered so that every graph belongs to (precisely) one vertex structure, 
and such that 



• the graphical details that are not needed for further computations are hidden, 



• equivalence classes £1X6 £IS large as possible, that is, combine as many graphs 
as possible. The number of vertex structures then becomes rather small. 

Summmation over graphs then is replaced by summation over vertex structures, 
with appropriate weights. We intent to make the set of vertex structures an algebra 
with operations that have their analog for graphs (such as concatenation), but oth- 
erwise are much more efficient and fast. In this algebra the series of tadpoles and 
the correlation functions are constructed. 

Once we know this structure it is straightforward to implement the construction 
of the series expansions in a language like C or C++. 

In Sect. 2 we recapitulate the principles of the linked cluster expansion in graph- 
ical terms. This is done only to the extent required to understand the necessary 
modifications for non-zero vacuum expectation values. We will show that this is 
conveniently implemented by a second kind of vertex renormalization, supplemented 
by a recurrence relation for computing the tadpole moments. In order to obtain 
comparable efficiency and order as in the symmetric case, in Sect. 3 the algebra of 
weighted vertex structures is introduced. The hopping parameter series construc- 
tion then is done in this algebra, generalizing the graphical expansion to a vertex 
structure expansion. A summary is given in Sect. 4. 

2 The graphical approach to series construction. 
2.1 Principles 

We start with a short survey of the series construction. We put particularly em- 
phasis on those points important to understand the differences of models with non- 
vanishing vacuum expectation value. For many details that are the same as for 
symmetric models and that are not repeated here, the reader is referred to [EJ. 

Let A denote a D-dimensional hypercubic lattice. It is assumed to be of shape 
Lq x L\ x • • • x Ld-i, with the Li > 4 even, finite or infinite, and periodic boundary 
conditons imposed in the directions of finite length. We discuss models described 
by a partition function of the form 



exp {W{J, v)) = Z(J, v)= \\ dfj,(p{x)) exp {-S hop {p, v) + ^ J{x)p(x)), (4) 



and 
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where p denotes a real valued scalar field, and J are external sources used to generate 
the correlation functions. The action Shop is assumed to be of the form 



S{p,v) = -- ^2 p(x)v{x,y)p{y). (5) 



2 

For pure nearest neighbour interaction, 

/ \ f 2k , x = y ± v for some u, 
V ^ y) = { , otherwise, (6) 

with v the unit vector in the positive vth direction, v — 0, . . . , D — 1. We confine 
attention to this case, although more general pair interactions are feasible M. The 
single site measure dp(p) is appropriately bounded for the partition function to 
exist and to be analytic for small hopping interaction v(x,y) and small external 
sources J(x). For instance, let I be some real interval, finite or infinite, and xi the 
characteristic function of /, 

{ \ otherwise. ^ 

Then the hopping parameter series discussed below are convergent for 

dp(p) = dp X i(p)e- v ^, (8) 

with V(p) bounded from below by 

V(p) >cp 2 -d, c> 0, (9) 

cf. [D, |J. The two examples given in the introduction correspond to the single site 
measures 

dp(p) = dp exp (-p 2 - X(p 2 - l) 2 + Hp) (10) 

and 

dfi(p) =dp p X[o,oo[(p) ex P (-V(p)) (11) 

(where p = £ 2 , cf. (0)), with V(p) bounded according to (Bp. Connected correlation 
functions are obtained by differentiation of the generating functional, 

W^(x 1 ,...,x n \v)=<p(x 1 )---p(x n ) > c = — — ——W(J,v) . (12) 

dJ(xx)---dJ(x n ) J=0 

Susceptibilities are defined as zero momentum correlations such as 

D-l 

xeA xeA i=o 
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Renormalized coupling constants are derived from them in a straightforward way. 

The linked cluster expansion is the Taylor expansion of the generating functional 
W(J,v) with respect to v(x,y), 



W{J,v) =(exp^ V ( X ,y)-JL- ) W (J,v) 



(14) 

v=0 



It generalizes in the obvious way to connected correlation functions. Multiple deriva- 
tives of W with respect to v(x, y) are managed by the identity |J 



dW 1 / d 2 W dW dW \ 

dv(x,y) ~ 2 \dJ(x)dJ(y) + dJ(x) dJ(y) ) ' 



(15) 



Correlation functions become series in the hopping parameter k, with coefficients 
depending on the toplogy of the lattice and the interaction /i, Eqn.(H). 

The evaluation of the Taylor expansion rapidly becomes complex with increasing 
order, which is the number of derivatives with respect to v(x, y). To manage it one 
is led in a natural way to a graphical device. In the present context, a graph is a 
strucutre 

T = (C r ,B r ,E r ^ r ), (16) 

where Cr and Br ^ are disjoint sets, the set of internal lines and the set of vertices 
of T, respectively. Er is a map 

E r :B r ^ {0,1,2,...}, u -> E r (u), (17) 

that assigns to every vertex u G Br the number of external lines Er(u) attached 
to it. If Er(u) 7^ 0, u is called an external vertex of T. The number of external 
lines of T is given by ^2 u&Br Er(u). Finally, $r is the incidence relation that assigns 
internal lines to their two endpoint vertices. Selffines are excluded. We consider 
lines as being unoriented, so $r maps onto unoriented pairs of vertices 



$r : £r -> (Br x B T ). (18) 



where for v,w G Br we have (v,w) = (w,v). For every pair v,w G Br, w/w, the 
number of common lines of v and w is defined by 



m(v, w) = $ r x ((f , w)). (19) 
The number of lines attached to u G Br is given by 

l(u) := m ( u > w ) + Er{u). (20) 

Connectedness is as usual path connectedness. 
Two graphs 

r 1 = {£ 1 ,B 1 ,E 1 ,<S> 1 ), T 2 = (£ 2 ,B 2 ,E 2 ,$ 2 ) (21) 
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are called (topologically) equivalent if there are two maps fa : B± — > B 2 and 2 : 
L\ — > £ 2 , such that 

$ 2 o0 2 = 0io$i, E 2 o ( j) 1 = E u (22) 
where o means decomposition of maps, and 



c/> 1 :B 1 xB 1 ^B 2 x B 2 , <t> x {v,w) = (fW^iH). (23) 

A symmetry of a graph T = (£, £>, E 1 , $) is a pair of maps <pi : £> — > B and </> 2 : £ — > 
£, such that 

$002 = ^0$, Eo(p 1 = E. (24) 

The number of those maps is called the symmetry number of T. 

The set of equivalence classes of connected graphs with E external and L internal 
lines is henceforth denoted by Q E {L), and 

Qe ■= |J Qe(L). (25) 

L>0 

Further graph classes will be introduced when they are needed below. Graphical 
expansions are done with respect to equivalence classes of graphs. In the following 
a graph T is understood to represent the complete class of graphs equivalent to V in 
the sense above. 

Connected correlation functions get a representation as a sum over connected 
graphs, each graph being endowed with the appropriate weight. In this descrip- 
tion a line of a graph represents a hopping propagator. For instance, the 2-point 
susceptibility becomes 

X2 = ^(^f E fl ( r ), (26) 
l>o reg 2 (L) 

a(T) = w(T) J] v l{u) fa), (27) 

u&B v 

where 

v i (/ J ,) = W {i \0,...,0\v = 0). (28) 
w(T) is a rational number that accounts for the various topological symmetry and 

o c 

lattice imbedding numbers of T. The couplings v i (fj,) are real numbers and depend 
only on the local coupling constants of the model, that is on the single link measure 
ix. They are obtained by numerical integration. 

The computation of correlation functions proceeds by constructing more re- 
stricted and hence smaller graph classes. They are then composed of the latter, ei- 
ther analytically or by concatenation. The first step is the introduction of 1-particle 
irreducible (1PI) graphs. A graph is called 1PI if it has the following property: 
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Remove a single line. The resulting graph then has only one connected compo- 
nent which has external lines attached. For instance, the 2-point susceptibility X2 
becomes 

xl PI 

X2 = i-( 2K ) 2 (2D) X r' (29) 

in graphical terms 




+ 



(30) 



with 

x\ pi = y,^) l E ( 31 ) 

l>o rgg 2 1FI (L) 

that is, x\ FI is composed of IPI graphs only. 

The second necessary step is to compose IPI graphs of the even more restricted 
classes of 1-vertex irreducible (1VI) graphs and renormalized moments. This is 
called vertex renormalization . A graph V is called 1 VI if the following condition 
is satisfied. Remove an arbitrary vertex u of T, together with the external and 
internal lines attached to u. We denote the resulting graph by r u . Then, every 
connected component of T u has at least one external line left attached to one of its 
vertices. We write 

S k (L) = {r G Gf\L) I r is 1VI } (32) 

for the set of graphs that are both IPI and 1VI. On the other hand, renormalized 
moment diagrams are IPI graphs that have exactly one external vertex, 



Qk{L) 



r has 1 and only 1 external vertex}. 



(33) 



The reason for the ~ will become clear in the next section. With both of these 
notions, susceptibilities are represented as 



T&S E (L) 



(34) 



L>0 



where the weights a(T) are computed in the same way as before by (PTI), with the 

oc 

following exception only. The vertex couplings v i (/i) therein are replaced by the 
renormalized moments 



(35) 



L>0 



1VI and renormalized moments provide a kind of product representation of IPI 
graphs by 1VI diagrams and renormalized moments. In the final end, after having 



S 



determined all the weights of the graphs, the coefficients are reorganized to give the 
series of xW 1 ■ 

Finally, we need to construct all of S k and They are generated from a 
sequence of appropriate base classes. Also, in order to do the graph generation 
on the computer, we need an algebraic representation of graphs. This is primarily 
done by so-called incidence matrices. The representation needs to be one-to-one 
and efficient, that is, fast algorithms are required to identify graphs by putting 
them into "canonical" representations. These problems are solved such that for 
O(N) symmetric models the most important series are available to 20th order in the 
hopping parameter. 

There are fast algorithms to generate 1-line irreducible (1LI) graphs. A graph is 
called 1LI if it cannot be divided in two by cutting a single line. For models with 
vanishing vacuum expectation value < p >= 0, a graph is 1LI if and only if it is 
1PI, so the notions 1LI and 1PI are equivalent for them. This concerns the graph 
classes S k and Q k . 

However, if the expectation value of the field < p >^ 0, the properties of being 
1LI and being 1PI are no more the same. 1LI implies 1PI, but not vice versa. The 
classes of renormalized moments Q k become enormously enhanced. This problem is 
addressed in the next subsection. 

2.2 Tadpoles 

Now suppose that the model under consideration is not subject to a Z2 symmetry. 
We consider the necessary modifications of the the series construction. 

• Correlation functions with an odd number of fields do no more vanish in gen- 
eral. This implies that higher connected correlation functions have a more 
complicated representation in terms of 1PI ones. 

• Graphs are no more restricted to having vertices only with an even number of 
lines attached. " Odd" vertices are allowed now. 

• The vacuum expectation value of the field does no more vanish in general, 

Xi =< p(x) > ^ 0. (36) 

The first two generalizations are easily included. First, the map between con- 
nected and 1PI correlations is done in closed analytic form. We are still mainly 
concerned with 1PI graphs only. Furthermore, vertices with odd number of lines do 
not provide a major problem. The inclusion of 1LI graphs having vertices with an 
odd number of lines attached implies an enormous enhancement of their numbers. 
For instance, the graph class S2(L) with number of internal lines L < 18 grows 
by about 3 orders of magnitude, as can be seen from Table 1 below. However, the 
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algorithms are so efficient that they still work with this generalization. They require 
a minimal modification only, as long as 1LI graphs are concerned. 

A non-vanishing one-point function fl36| ) provides the major problem. Graphi- 
cally we represent this expectation value by the "tadpole", 

Xi = < p{x) > = {~) 

(37) 

Now 1PI graphs are no more necessarily ILL In particular, by the very definition, 
tadpole graphs are always 1PI, but almost never ILL To obtain the graphical rep- 
resentation of 1PI n-point functions for n > 2 we have to 

• find a way to construct the tadpole graphs, and 

• attach the tadpoles to 1LI graphs in all possible ways. 




Z-^m>0 ml 




(38) 

We can do better. In the last subsection we have recapitulated the composition 
of 1PI graphs by 1VI graphs and renormalized moments. Whereas the 1VI graphs 
of Sk are already 1LI, as follows from their definiton, the renormalized moments 
Qk are 1PI only. Hence it is sufficient to do the above tadpole attachment only to 
generate Qk- This amounts to a further vertex renormalization. We circumvent the 
construction of 1LI graphs with large number of external lines, which would be very 
expensive. 

We define the sets of first renormalized moment graphs Qk(L) by 

Q k (L) = {T E Gl L1 (L) | L has 1 and only 1 external vertex}, (39) 
and the first renormalized moments themselves by 

XQk = £(2/s) L E w ^ II "?(«) (A*)- (40) 

L>0 reQfe(L) ueB r 

The first renormalized moment graphs are 1LI and are constructed as for Z 2 sym- 
metric models. The second renormalized moments diagrams Qk, O), are then 
obtained by attaching tadpoles to the vertices of the graphs of Qk- Let us assume 
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for the moment that the series representation of the tadpole moment Xi is known, 
for instance in the form 



= J2( 2k ) L Xi(L). (41) 



L>0 

For convenience we introduce the quantities 

Xi («,//) = ^2v(x,0) Xi{K,fj) = (2k) (2D)xi(k,h), 

xfW) = -^(xfW))" = (2k)" ^l( Xl ( K ,^)y, u>2. (42) 
With the tadpole vertices v^ ta,dp (K, \i) defined by 

■ c ' tadp f K , n) =vl fa) + £ v° i+v (/,) m), (43) 



the second vertex renormalization (^) is replaced by 

£(a*)-'*(k,/0 = £(2k) £ £ «(F), (44) 
i>o reQ,(L) 

with 

5(r) = «;(r)n<(5 dp (^^)- (45) 

Finally, after having determined all the weights of the graphs, the coefficients are 
reorganized as before to give the series of xW 1 ■ 

We have assumed that the hopping parameter series ( |4~l]) of the 1-point function 
is known. Actually, this series representation is obtained by (|44j) itself for i — 1. 

For, we have i>i(k, /i) = xi( K ? / i ) ? an d Xi"^?/*) = 0(k 1/ ). Hence (f42|-[45|) becomes 
the recursion relation 

l>o reQi(L) «eB r V v>i / 

(46) 

o c 

with initial condition xi(^,/-0 =v 1 (/i) + O(k). This recursion is solved order by 
order in k. We need the first renormalized moments Qi(L) with one external line, 
consisting of 1LI graphs only 

In summary, the generation of the hopping parameter series of m-point functions 
to order M, say, according to the discussion above proceeds as follows. 

1. The single site coupling constants v i (/x) for 1 < i < m + M are computed. 
This amounts to the evaluation of a sequence of low- dimensional integrals. 
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2. Construction of the graphs of Se{L) for E <= m and L < M, and of the 
first renormalized moments Qe{L) for E + L < m + M, alltogether being 1LI, 
and their weight factors. This is done as for Z2 symmetric models. The only 
generalization here is that we allow for vertices having an odd number of lines 
attached. This is easily included. 

3. The series of the tadpole moment xi are computed by the recursion relation 
( |4q) , starting with xi\ K i A 4 ) =v i (fv + 0(k). This requires the knowledge of 
the first renormalized moments Qi(L) for L < M, obtained in step 2. 

4. Once the series of Xi are known to order M, we compute the second renormal- 
ized moments v°(k, fj) according to ( [i"3"[{4li)| ) , by attaching the tadpole moment 
in all possible ways to the vertices of the first renormalized moment graphs 
(the first vertex renormalization). These second renormalized moments then 
are attached to the vertices of the the graphs of Se{L) with E < m and 
L < M (the second vertex renormalization). In this way we obtains the series 
coefficients of the desired 1PI correlation functions. The final computation of 
the connected correlations is done analytically. Of course each of the steps 
above involves a sorting of the series coefficients according to their order in k. 

3 Vertex structures. 
3.1 Motivation 

As we have seen in the last section, the only graphs that are actually generated 
and operated on are always ILL Their numbers for the simplest classes are listed 
in Table 1. By the procedure outlined at the end of the last section, all the other 
moments and correlations are no more obtained as a graphical representation but 
directly as hopping parameter series by the two vertex renormalizations. 

Although this construction works in principle, there is need for improvement. 
This is for the following reasons. 

• Beyond exceptional cases, in most applications the coupling constants are not 
fixed numbers but vary over some subsets of the reals, for instance for the 
determination of critical surfaces and universality domains. For every point 
in coupling constant space, the series construction to 18th order as described 
above will take a couple of hours on a Sparc Sun work station. For many 
points this becomes too expensive. 

• The 1LI graphs generated are needed for each computation and hence have to 
be stored. Even in binary coded form this needs a couple of GByte. 

• The most severe point is the ill-conditionedness of the composition procedure. 
If the weights of the graphs involved would be rational (as e.g. for the Ising 
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Table 1 : The number of inequivalent graphs for the various 1-line irreducible graph 
classes defined in Section 2. For the first renormalized moments Qe{L) this number 
does not depend on E, the number of external lines. For comparison we have given 
also the number of graphs having only vertices with even number of lines attached. 
They are denoted by the superscript cv . Notice that 5| v = 0. 



L 


Qf{L) 


Qi{L) 


S?(L) 


S 2 (L) 


Sa(L) 





1 


1 


1 


1 


1 


1 
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o 


o 


o 


o 


2 


1 


1 





1 


1 


3 





1 


1 


1 


1 


4 


4 


4 





4 


6 


5 





6 


2 


6 


12 


6 


15 


23 


3 


23 


52 


7 





49 


8 


53 


145 


8 


79 


174 


9 


190 


567 


9 





483 


40 


575 


1954 


10 


439 


1681 


68 


2089 


7577 


11 





5446 


247 


7349 


28811 


12 


2877 


19562 


470 


27908 


115290 


13 





69575 


1779 


106666 


465112 


14 


20507 


260752 


3937 


423687 


1930270 


15 





989445 


14801 


1710920 


8134945 


16 


161459 


3874201 


35509 


7080691 


34982052 


17 





15446928 


135988 


29828057 


152884073 


18 


1376794 


63028306 


350614 


128100002 





model in zero external field), this problem is absent. However, whereas the 

°° 

weights w(T) are rational numbers, in most cases the vertex couplings v i (/i) 
are non-rational. They are of both signs, of comparable size, and they are 
computable to limited precision only. We expect that roundoff errors rapidly 
accumulate, limiting the series computation above to about the 10th order in 
the hopping parameter. 

A solution to this problem is provided by expansion into vertex structures. Ver- 
tex structures have already been introduced and used for Z 2 symmetric models pi . 
They provide the possibility to arrange the rational parts of the weights according 
to the precsribed concatenation and recursion relations, keeping track of where and 
how the vertex couplings v i (/x) are to be multiplied and added. The actual mul- 
tiplication with the vertices is done only in the very end, when all the other terms 
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have been collected, to obtain the final series coefficients of the correlation functions. 
At this stage then this amounts to a small number of additions and multiplications 
(a couple of thousands). Hence, using operations on vertex structures, the series 
construction will become well-conditioned even to large orders. Also, instead of 
storing the graphs of <S& and Qk, one only needs to keep the vertex structures or 
their algebraic representations, which are both tiny objects and small in number. 

3.2 Definitions and tadpole moment construction 

In order to work with vertex structures and to understand in which sense they 
replace graphical manipulations we need to introduce some basic definitions. 
A vertex structure v is an ordered sequence of non-negative integers 

v = (z^)igN = 0i,^2,---); ^ G {0, 1,2, . . .}, (47) 

with only finitely many z/j ^ 0. For v given by fl47|) we call Vi = v^ the ith component 
of v, s(v) = max„^o % the size or the length of v, and n(v) = J2i>i v i i the number 
of lines of v. The reason for this becomes clear below where a vertex structure is 
associated to a graph. We always have s(v) < n(v). The set of all vertex structures 
v with n(v) = n is denoted by V n . We put 

V (n) = |J V fc . (48) 

l<fc<n 

For non-negative integer k we define the vertex structure k G Vk by 

k = (fi)ieN) with v k = 1 and v^ = for % ^ k. (49) 

To a given graph Y a vertex structure v = v(T) is associated as follows. We set 

v(V) = (vx,v 2 ,...), 

Vi = \{u E B r | l{u) = i}\, (50) 

that is, Vi is the number of vertices of T that have precisely i lines attached. If F has 
E external and L internal lines, we have n(v(T)) = 2L + E, that is, v (T) G V<il+e- 

For given T, the vertex structure v (T) is uniquely defined, but not vice versa. 
For given vertex structure v there are many graphs T such that v = v(T). Their 
number is large but finite. This allows us to define a very useful equivalence relation 
~V5 in the graph classes under consideration. Hereby, two graphs Ti and T 2 are 
equivalent, Ti T 2 , if v(Ti) = v(T 2 ). This equivalence relation devides the graph 
classes into small sets of equivalence classes, each class being represented by a single 
vertex structure. We will associate appropriate rational weights to vertex structures 
and will define operations on them that are equivalent to the concatenation rules 
and recursion relations of the last section. 
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Toward this end we endow U n >iV ( - Tl - ) with the structure of an associative algebra. 
First, for later convenience we define the product of two vertex structures v and v' 
by 

v-v' = (Vi + u/) i6 N- (51) 

By some "external" process, such as the graph association above, a vertex struc- 
ture gets associated a weight. A pure, weighted vertex structure is a pair (w v ,v) 
of a number w v of some field JF, for our purpose the field of rational numbers, 
and a vertex structure v. For w v = 1 we use the notation V v = (l,v) and write 
w v • V v = (w v , v). A weighted vertex structure then is a finite sum of the form 

(w v ) veV (n) = ^ w vV v . (52) 

This set can be given the structure of an T vector space in the obvious way, denoted 
by W^ n \ Furthermore, for 1 < k < n become sub- vector spaces of WW, which 
we denote by Wk- Finally, with Wo = J 7 , 

W = 0W„ (53) 

n>0 

is made a graded algebra with unity by defining for v G W s , v' G Wt the product 
V v • V V ' = P„.„' G W s+t , and by linearity. 

It is now easy to translate the series construction into this algebra. We examplify 
this for the computation of the tadpole moment. First, for the simpler graph classes 
(the 1LI ones), summation over graphs is replaced by summation over associated 
vertex structures. This is done by adding the rational weights of all graphs that 
belong to the same equivalence class of ~vs- For instance, let us define for L = 
0,1,2,... 

Qi{L) = {v G V 2L+1 | v(T) = v for some T G Qi(L)}. (54) 
For the first renormalized moments with series representation 

XCi = ^(2^) L giW, (55) 

L>0 

the coefficient of (2k) l is now obtained as 

s(v) 

qi (L)= ™(n n 5 «(«)^)= e MQi)n(K(ri) v \ (56) 

reQi(L) ueB v veQi(L) i=1 

where for every vertex structure v G Qi(L) 

TeQi(L)withv(T)=v 
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The latter sum consists of a large number of terms, but the w(T) are rational num- 
bers. On the contrary, in fl56|) the various terms are real numbers only, but their 
number is very small. With the identification of 



n 

W ■ 

i=l 



n 

n ifi (^r - w ■ v ^ v = ^ ^ ■ ■ o» (58) 



qi(L) becomes replaced by 



q x (L) = w v(Qi)-^ e W. (59) 

veQi{L) 

D-l 



To the lowest orders, we have e.g. on the 4 x oo lattice 



9i (0) 


= i-^w, 


ft(l) 


= o, 


ft (2) 


= ^-^(0,1,1), 


ft (3) 


— J ' ' (0,0,1,1) ; 


ft (4) 


= ^ ■ ^(o,o,o,i,i) + (QD 2 
D 2 

+ ' H0,2,0,0,l)> 



(60) 



where trailing zeros of the vertex structures have been cut off. The recursion for the 
series construction of the tadpole moment Xi{ K -> aOj Eqn. PB]), is converted into the 
recursion in the algebra W as follows. Let us set 

Xi(L = 0) = 7> T , (61) 

which is the translation of the initial condition of fl46|). The analog of (H^) becomes 

X?{L) = 2D-xi(i-l), ^>1, 



X?(L) = -■ jr ^\L').^r l \L-L'), 2<u<L. (62) 



v 

L'=l 

The tadpole moment recursion in W then reads for L = 1,2,... 

ul) = E E <Qi)-nft E (^.o-n+E^-x^^; 

i'^ueQi^') i=lj=l L,j=0 y i/=l 



(63) 
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where the prime indicates that the summations are restricted to J^-. Ly — L — L' . 
The solution of this recurrence relation yields up to some order M, say, the vertex 
structure representation of Xx(L), that is, 

Xi(L) = w v(Xi)-Pv, 0<L<M. (64) 

To the same order as for qi above we obtain, again on the 4 x oc D ~ l lattice, 

Xl(0) = 1 •?>(!), 

Xi(l) = 2D-V ( x,x), 

Xi(2) = D-V i0 ,i,i ) +4D 2 -V {1>2) + 2D 2 -V m i), 

Xl(3) = y • 7>,0,1,1) + 21)2 • ^(1,0,2) + ■ ^(1,1,0,1) + 2# 2 • ^(0,2,1) 

4/13 

+ 8D 3 ■ V m + 12D 3 - P (2jl)1) + — ■ 7^(3,0,0,1) , (65) 

Xi(4) = (6D 2 - 3D + 1) • P (0j3i1) + — • P (0 , 2 ,o,o,i) + ^ ■ ^(o,o,o,i,i) 
2D 2 2D 2 , „ 

+ -y • ^(1,0,0,2) + -y • ^(1,0,1,0,1) + QD ■ ^(2,0,1,1) + 12 ^ 3 • V(l,l,2) 

5D 2 

+ SD 3 • P(i, 2 ,o,i) + 2D 3 • P (2 ,i, 0> o,i) + -y ■ ^(0,1,1,1) + 16Z> 4 • P(i,4) 

4 32£> 4 4 2L> 4 

+ 48D • V(2,2,i) H g — • V(3 t x ) o 1 x) + SD ■ "P(3,o,2) H — g— • ^(4,0,0,0, i)- 

The final series representation of the tadpole moment Xi? (PI) , is obtained by the 
map fl58f) , that is, 



s (") !J . 

L>0 ^eV 2 L+i »=i 



In a similar way as the computation of the tadpole moments, the vertex renormal- 
ization processes are done in terms of vertex structures, that is in the algebra W. 
This provides a well conditioned algorithm. To order 18 we have about 5 orders 
of magnitude less numbers of vertex structures to deal with than we would have 
graphs. In Table 2 we summarize some numbers on vertex structures. 



4 Summary 

We have generalized the linked cluster expansion to models with nonvanishing vac- 
uum expectation value. The main emphasis here was on the technical details neces- 
sary to achieve the hopping parameter series of correlation functions to comparable 
order as for symmetric models. The major additional problems that had to be solved 
are the following ones. 
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Table 2: The number of mutually different vertex structures associated with various 
graph classes according to used for the computation of the hopping parameter 
expansion to order 18. Compared to the number of graphs as given in Table 1, 
there are about 5 orders of magnitude less number of vertex structures than graphs. 
TZi(L) denotes the vertex structures of the tadpole moment \i t° order L, that is 
of the 1-point function < p(x) >. 



L 


Qi(L) 


S 2 (L) 


S 3 {L) 


fti(L) 


o 


1 


1 


1 


1 


1 


n 


n 


n 


1 


2 


1 


1 


1 


3 


3 


1 


1 


1 


7 


4 


4 


3 


5 


15 


5 


6 


5 


8 


31 


6 


15 


12 


19 


60 


7 


22 


19 


30 


110 


8 


45 


37 


55 


198 


9 


69 


58 


89 


342 


10 


116 


103 


148 


569 


11 


183 


161 


229 


943 


12 


294 


262 


364 


1513 


13 


431 


399 


551 


2377 


14 


677 


628 


839 


3700 


15 


1003 


938 


1246 


5672 


16 


1470 


1413 


1838 


8529 


17 


2148 


2068 


2676 


12745 


18 


3119 


3047 




18775 



• The tadpole moment representing the series expansion of the vacuum expec- 
tation value < p(x) > has to be worked out. It obeys a nonlinear recursion 
relation that involves 1-line irreducible moments only. 

• Computation of connected correlation functions requires the attachment of 
tadpole moments to all vertices. This is conveniently done by including a 
second vertex renormalization. 

• There is a huge proliferation of graphs and in turn of roundoff errors due to the 
ill-conditioned real arithmetic for generic coupling constants. At the first sight 
this limits the computations of the series to order 10 in the hopping parameter. 
The solution of this problem is to replace the graphical expansion of correlation 
functions and moments by an expansion in terms of weighted vertex structures. 
The latter can be given the structure of a graded algebra W with unity over the 
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rational numbers. The hopping parameter series construction is done mainly 
in this algebra. 

The vertex structures and their weights have to be computed only once. For 
the most important moments they are available to 18th order. Using them, for any 
set of coupling constants (the single site measure (§)), the hopping parameter series 
themselves are obtained within a couple of minutes on a Sparc Sun workstation. 

We have dicussed real valued fields only, merely for simplicity. The ideas behind 
vertex structures are generalized in a straightforward way to models with larger 
symmetry group such as O(N), broken down e.g. to 0(N — 1). This only requires 
slight changes of equivalence relations and extensions of the algebra W. 

Various applications are in preparation and will be presented elsewhere. 
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